Climate Forecast System Version Research Articles

Midsummer precipitation prediction over eastern China by the dynamic downscaling method

AbstractThis study assesses the midsummer precipitation prediction over eastern China by the dynamic downscaling method. Based on the Climate Forecast System version 2 of the National Centers for Environmental Prediction and the Weather Research and Forecasting Model, the prediction performance of global and regional models on the July precipitation over eastern China is further analyzed by hindcast experiments from 1982 to 2010 and prediction experiments from 2011 to 2021. The results suggest that the regional model forced by the global model can noticeably improve the prediction skill for precipitation in eastern China, especially in the region from the South of North China to the Yangtze River Basin, referred as the northern China in this paper. In addition, we perform a diagnostic analysis of the reason for the improvement of the model prediction skill. The results indicate that the high resolution of the regional model and the refinement of physical process parameterizations contribute to improving the simulation ability for the East Asian atmospheric circulation pattern, heat flux, especially for the meridional teleconnection pattern in East Asia and the sensible heat flux in the northern China, thus further improving precipitation prediction.

Synoptic Analysis and Subseasonal Predictability of an Early Heatwave in the Eastern Mediterranean

Greece and the surrounding areas experienced an early warm spell with characteristics of a typical summer Mediterranean heatwave in mid-May 2020. The maximum 2 m temperature at Kalamata (southern Greece) reached 40 °C on May 16th and at Aydin (Turkey), it was 42.6 °C on May 17th. There was a 10-standard deviation positive temperature anomaly (relative to the 1975–2005 climatology) at 850 hPa, with a southwesterly flow and warm advection over Greece and western Turkey from May 11 to 20. At 500 hPa, a ridge was located over the Eastern Mediterranean, resulting in subsidence. The aims of this study were (a) to investigate the prevailing synoptic conditions during this event in order to document its occurrence and (b) to assess whether this out-of-season heatwave was predictable on subseasonal timescales. The subseasonal predictability is not a well-researched scientific topic in the Eastern Mediterranean Sea. The ensemble global forecasts from six international meteorological centres (European Centre for Medium-Range Weather Forecasts—ECMWF, United Kingdom Met Office—UKMO, China Meteorological Administration—CMA, Korea Meteorological Administration—KMA, National Centers for Environmental Prediction—NCEP and Hydrometeorological Centre of Russia—HMCR) and limited area forecasts using the Weather Research and Forecasting model with the Advanced Research dynamic solver (WRF) forced by Climate Forecast System version 2 (CFSv.2; NCEP) forecasts were evaluated for lead times ranging from two to six weeks using statistical scores. WRF was integrated using two telescoping nests covering Europe, the Mediterranean basin and large part of the Atlantic Ocean, with a grid spacing of 25 km, and Greece–western Turkey at 5 km. The results showed that there were some accurate forecasts initiated two weeks before the event’s onset. There was no systematic benefit from the increase of the WRF model’s resolution from 25 km to 5 km for forecasting the 850 hPa temperature, but regarding the prediction of maximum air temperature near the surface, the high resolution (5 km) nest of WRF produced a marginally better performance than the coarser resolution domain (25 km).

A hybrid statistical downscaling prediction model with timescale decomposition for summer precipitation over Hainan Island, China

AbstractBy observing well‐known sea surface temperature (SST) indices and outputs from Climate Forecast System version 2 (CFSv2) hindcasts and predictions, in this study, we developed a hybrid statistical downscaling prediction model (HSDPM) based on a timescale decomposition approach, which effectively improved the forecasting ability of summer (June–August) precipitation over Hainan Island (HNSP) in China; this precipitation presented multiple obvious timescale variabilities. We found that the interannual variability of the HNSP is closely related to the observed North Atlantic tripole SST anomalies in the preceding winter and the summer sea level pressure over the South China Sea by CFSv2 predicted in March. On the other hand, the interdecadal variability in the HNSP is linked to the prewinter Pacific Interdecadal Oscillation (PDO) according to the observations, the CFSv2 predictions of the summer surface air temperature over the tropical central‐east Pacific. On this basis, both the interannual and interdecadal components of the HNSP are effectively predicted using the corresponding predictors via multiple regression; thus, the HSDPM is ultimately established by combining the above two components. Compared with the original CFSv2 model, the HSDPM model achieves a considerable improvement in performance in predicting the HNSP. Specifically, the temporal correlation coefficient, ratio of root‐mean‐square error and anomaly sign consistency rate between the observed and HSDPM‐predicted HNSPs are 0.70, 17.4% and 80.5%, respectively, which are significantly greater than the above three forecasting index values of 0.23, 38.3% and 48.8%, respectively, obtained from the original CFSv2 predictions. The application of the HSDPM may be beneficial for drought and flood prevention and mitigation in Hainan Island.

Effects of Arctic sea-ice concentration on turbulent surface fluxes in four atmospheric reanalyses

Abstract. A prerequisite for understanding the local, regional, and hemispherical impacts of Arctic sea-ice decline on the atmosphere is to quantify the effects of sea-ice concentration (SIC) on the turbulent surface fluxes of sensible and latent heat in the Arctic. We analyse these effects utilising four global atmospheric reanalyses, ERA5, JRA-55, MERRA-2, and NCEP/CFSR (including both the NCEP Climate Forecast System Reanalysis (CFSR) and the NCEP Climate Forecast System Version 2 (CFSv2)), and evaluate their uncertainties arising from inter-reanalysis differences in SIC and in the sensitivity of the turbulent surface fluxes to SIC. The magnitude of the differences in SIC is up to 0.15 but typically around 0.05 in most of the Arctic over all four seasons. Orthogonal-distance regression and ordinary-least-squares regression analyses indicate that the greatest sensitivity of both the latent and the sensible heat flux to SIC occurs in the cold season, November to April. For these months, using daily means of data, the average sensitivity is 400 W m−2 for the latent heat flux and over 800 W m−2 for the sensible heat flux per unit of SIC (change in SIC from 0 to 1), with differences between reanalyses that are as large as 300 W m−2 for the latent heat flux and 600 W m−2 for the sensible heat flux per unit of SIC. The sensitivity is highest for the NCEP/CFSR reanalysis. Comparing the periods 1980–2000 and 2001–2021, we find that the effect of SIC on turbulent surface fluxes has weakened owing to the increasing surface temperature of sea ice and sea-ice decline. The results also indicate signs of a decadal-scale improvement in the mutual agreement between reanalyses. The effect of SIC on turbulent surface fluxes arises mostly via the effect of SIC on atmosphere–surface differences in temperature and specific humidity, whereas the effect of SIC on wind speed (via surface roughness and atmospheric-boundary-layer stratification) partly cancels out in the turbulent surface fluxes, as the wind speed increases the magnitudes of both upward and downward fluxes.

Unusual subseasonal variability of Indian summer monsoon rainfall in 2020

AbstractAbove‐normal all‐India summer monsoon (ISM) rainfall was received in 2020 with 109% of its long‐period average. According to the India Meteorological Department (IMD), ISM 2020 experienced robust month‐to‐month rainfall variations; in particular, the peak monsoon month of July 2020 witnessed deficit rainfall. The present study investigated the possible mechanisms that are responsible for the July 2020 deficit rainfall. Analysis revealed that the strong low‐level moisture divergence caused by the westward‐propagating atmospheric Rossby wave, induced by suppressed western North Pacific (WNP) convection, is primarily responsible for the July rainfall deficit. The WNP suppressed convection is linked to anomalous low‐level anticyclonic circulation. It is noted that the intense WNP anticyclone is maintained by the tropical Indian Ocean (TIO) warming‐induced atmospheric Kelvin waves and strong low‐level convergence over the Meiyu‐Baiu rainband. Unlike July, the strength of the WNP anticyclone and TIO warming are weaker in June and August. In addition to the observational analysis, the skill of the Climate Forecast System version‐2 (CFSv2) model in predicting the July rainfall at a lead of about two months is examined. Strong positive rainfall anomalies over India are predicted by the model in July, which is in contrast to the observations. Despite predicting the TIO warming, the model failed to capture the westward extension of the WNP anticyclone, which impacted the July rainfall prediction. To understand the importance of the WNP anticyclone on ISM rainfall, a CFSv2 sensitivity experiment is carried out by imposing strong negative sea surface temperature (SST) anomalies over the WNP region to excite the anticyclone. Both the westward extension of the WNP anticyclone and suppressed rainfall over the core monsoon region in July are captured by the model when low‐level circulation is imposed. The model sensitivity experiment strongly supports the role of the WNP anticyclone and its westward extension in inducing reduced rainfall in July 2020 over the monsoon trough region.

Data-Driven Global Subseasonal Forecast for Intraseasonal Oscillation Components

As a challenge in the construction of a “seamless forecast” system, improving the prediction skills of subseasonal forecasts is a key issue for meteorologists. In view of the evolution characteristics of numerical models and deep-learning models for subseasonal forecasts, as forecast times increase, the prediction skill for high-frequency components will decrease, as the lead time is already far beyond the predictability. Meanwhile, intraseasonal low-frequency components are essential to the change in general circulation on subseasonal timescales. In this paper, the Global Subseasonal Forecast Model (GSFM v1.0) first extracted the intraseasonal oscillation (ISO) components of atmospheric signals and used an improved deep-learning model (SE-ResNet) to train and predict the ISO components of geopotential height at 500 hPa (Z500) and temperature at 850 hPa (T850). The results show that the 10–30 day prediction performance of the SE-ResNet model is better than that of the model trained directly with original data. Compared with other models/methods, this model has a good ability to depict the subseasonal evolution of the ISO components of Z500 and T850. In particular, although the prediction results from the Climate Forecast System Version 2 have better performance through 10 days, the SE-ResNet model is substantially superior to CFSv2 through 10–30 days, especially in the middle and high latitudes. The SE-ResNet model also has a better effect in predicting planetary waves with wavenumbers of 3–8. Thus, the application of data-driven subseasonal forecasts of atmospheric ISO components may shed light on improving the skill of seasonal forecasts.

Observational characteristics and causes of extreme heavy rainfall in typhoon Duksuri (2305)

Super Typhoon Duksuri (2305), made landfall in Fujian Province, causing significant impacts from July 27th to July 31st. The typhoon resulted in heavy to torrential rainfall in multiple provinces, particularly in southern Taiwan, eastern Fujian, and the Beijing-Tianjin-Hebei region. Daily accumulated precipitation reached a maximum of 200 to 220 millimeters. Using ERA5 reanalysis data and NCEP Climate Forecast System Version 2 reanalysis data, this article conducted an observational analysis of Duksuri. The results indicate that the interaction between subtropical high-pressure systems, typhoons, and the strong southeast and east winds created favorable conditions for the prolonged duration of Duksuri and extreme precipitation. The extreme rainfall in southern Taiwan and eastern Fujian was primarily caused by the abundant moisture associated with the typhoon itself. In contrast, the extreme rainfall in the Beijing-Tianjin-Hebei region was primarily due to the northward movement of the remnants of the typhoon, hindered by the high-pressure systems, resulting in longer and concentrated precipitation. Additionally, both regions are located on the windward slopes of mountainous areas, which significantly enhance the rainfall due to orographic lifting.

Assessment of Extreme Ocean Winds within Intense Wintertime Windstorms over the North Pacific Using SMAP L-Band Radiometer Observations

Here, we examine extreme ocean winds associated with intense wintertime extratropical windstorms over the North Pacific. The study was mainly based on NASA Soil Moisture Active Passive (SMAP) L-band radiometer observations allowing the retrieval of ocean wind speeds up to 70 m/s regardless of precipitation intensity. Additionally, we assessed the ability of atmospheric reanalysis ERA5 and the Climate Forecast System Version 2 (CFSv2) to reproduce high-wind features within severe windstorms, particularly those associated with “explosive” cyclogenesis. The analysis identified 145 windstorm events with hurricane-force (HF) wind zones within the SMAP L-band radiometer swath from 2015 to 2023. These windstorms develop most frequently over two areas: southeast of Kamchatka and south of Alaska, spanning 40–47°N latitudes. Both reanalysis datasets significantly underestimated HF wind speeds compared to SMAP measurements, but CFSv2 tends to reproduce more-intense windstorms than ERA5. Among the notable new findings is that the SMAP data revealed two distinct groups in maximum wind speed distribution, indicating the existence of a separate class of severe windstorm events with a distinct mechanism for extreme wind formation related probably to a Shapiro–Keyser cyclogenesis and the presence of sting jet (SJ) feature. The study highlights the potential of SMAP measurements to study wind extremes and underscores the need for improvements in operational predictive models to better reproduce the formation of SJ windstorms.

Wave–Tide Interaction by Typhoon Ampil on Wave and Storm Surge in the Changjiang River Estuary and Its Adjacent Coastal Areas

The study used the SCHISM ocean model combined with the WWM III wind wave model to quantify the interaction between wind waves and tides in the coastal zone of the Changjiang River Estuary and its adjacent areas. The wave and storm surge during Typhoon Ampil, which made landfall in July 2018 in Shanghai, were simulated by using the Climate Forecast System Version 2 (CFSv2) and Medium-Range Weather Forecasts (ECMWF) latest reanalysis (ERA5) wind dataset from 1 July to 31 July. Model results with CFSv2 forcing show better performance in terms of significant wave height and storm surge than those with ERA5 forcing. To investigate the interactions between waves, water levels, currents, and their combined effects on significant wave and surge variations, six numerical sensitivity experiments were designed according to the different coupling methods between SCHISM and WWMIII. The research shows that in coastal areas with water depths of less than 10 m, waves are affected by water levels and currents. The differences in the effect on significant wave height between wave-tide coupling and one-way coupling with water levels and currents are negligible. Wave setup is an important physical term which cannot be ignored during the variations of storm surge caused by Typhoon Ampil. The contributions of wave set up were concentrated in coastal areas with water depths less than 10 m. The peak wave setup occurred in the Changjiang River Estuary, reaching 0.15 m. In Xiangshan Bay and Sanmen Bay, wave radiation stress makes the proportion of wave setup to the total surge reached more than 30%. The consideration of wave-tide interaction can effectively improve the accuracy of numerical wave and storm surge simulations, which can provide more accurate hindcasts of wave and storm surge variations in the Changjiang River Estuary and its adjacent coastal areas.

Predictability of Northwest Pacific marine heatwaves in summer based on NUIST-CFS1.0 hindcasts

Marine ecosystems in the Northwest Pacific (NWP) are being threatened by marine heatwaves (MHWs), which are abnormally high sea-surface temperature (SST) lasting for more than 5 days. However, there is limited research on the predictability of these MHWs. In this study, we investigate the predictability of NWP MHWs in summer using daily data from the coupled climate forecast system (version 1.0) of the Nanjing University of Information Science and Technology (NUIST-CFS1.0) for the period of 1984–2020. We demonstrate that the NUIST-CFS1.0 hindcasts can predict the spatial pattern of the total MHW days (HWD) over the NWP in summer with a lead time of up to eight months. Moreover, the linear trend and interannual variability can be predicted at lead times of up to nine and three months, respectively. We reveal that the interannual variability of HWD is strongly correlated with the preceding El Niño-Southern Oscillation (ENSO); and its predictability also shows significant interannual variability, with MHWs being more predictable in the decaying years of stronger El Niño events. In the El Niño decaying summers, there is an increase in the occurrence of MHWs; the spatial distribution of HWD can be predicted at lead times of up to 19 months, and its forecasting accuracy is far superior to that of the other years. The occurrence of MHWs during El Niño decaying summers is closely linked to strengthened NWP subtropical anticyclone (NWPSA), which increases solar radiation reaching the ocean surface and suppresses surface evaporation, leading to increased SST and thereby MHWs. The remarkable predictability of the MHWs in the NWP region during El Niño decaying summers can be attributed to the high skill at predicting ENSO and its associated NWPSA anomalies.

Predictability of the unprecedented 2022 late summer Yangtze River Valley and Tibetan Plateau heatwaves by the NCEP CFSv2

By using the one-month‑lead predictions from the Climate Forecast System version 2 (CFSv2), this study evaluates the predictability of the unprecedented heatwaves hitting the Yangtze River Valley (YRV) and Tibetan Plateau (TP) in 2022 late summer. CFSv2 successfully reproduces the central China heatwaves, but shows distinct performances between the YRV and TP. Specifically, CFSv2 exhibits stronger ability to capture TP heatwaves that is close to the highest historical prediction record, whereas the predicted extremity of the YRV heatwaves is underestimated compared to the observations. The accurate simulation of the intense upper-level circumglobal teleconnection (CGT) pattern, particularly for the East Asia CGT center (CGT EA), is essential for the success predictions of high temperature anomalies, due to the strong and accurate linkage between CGT EA and YRV/TP heatwaves in CFSv2. The intense CGT pattern is associated with a stronger model response to developing La Niña in CFSv2, quite different from the observed leading contribution from the unprecedented Pakistan rainfall. The impacts of developing La Niña and Pakistan rainfall on the intense CGT pattern are significantly overestimated and underestimated by CFSv2, respectively. Additionally, the La Niña-forced enhancement of lower-level western Pacific subtropical high (WPSH) plays another important role in YRV heatwaves, but the reproduction of WPSH intensity is unable to reflect YRV heatwaves in CFSv2 due to the southeastward shift of the predicted WPSH location, which should be responsible for the underestimated YRV heatwaves. The model improvements in aspect of the Pakistan rainfall, El Niño-Southern Oscillation–CGT linkage and lower-level WPSH location should help effectively predict the central China heatwaves in the future.

Intercomparison of tropospheric and stratospheric mesoscale kinetic energy resolved by the high‐resolution global reanalysis datasets

AbstractWith the development of advanced data assimilation and computing techniques, many modern global reanalysis datasets aim to resolve the atmospheric mesoscale spectrum. However, large uncertainties remain with respect to the representation of mesoscale motions in these reanalysis datasets, for which a clear understanding is lacking. The aforementioned challenges have served as a strong motivation to reveal and quantify their mesoscale differences. This study presents the first comprehensive global intercomparison of the tropospheric and stratospheric mesoscale kinetic energy and its spectra over two selected periods of summer and winter events among six leading high‐resolution atmospheric reanalysis products: European Centre for Medium‐Range Weather Forecasts (ECMWF) Reanalysis v5 (ERA5), China Meteorological Administration Reanalysis, Modern‐Era Retrospective Analysis for Research and Applications version 2, National Centers for Environmental Prediction's Climate Forecast System version 2 (CFSv2), Japanese 55‐year Reanalysis, and ECMWF Reanalysis‐Interim. A state‐of‐the‐art global operational model is adopted as a supplementary reference. Although all reanalysis datasets can reproduce broad distribution characteristics that are grossly consistent with the 9 km model, there are substantial discrepancies among them in magnitudes. The ability to capture mesoscale signals is closely linked to their resolutions, but it is also impacted by other factors, including, but not limited to, the selected types of energy, seasons, altitudes, latitudes, model diffusions, parametrization schemes, moist condition, assimilation methods, and observation inputs. Moreover, all datasets illustrate conclusive behaviors for the prevalence of the rotational component in the troposphere, whereas only very few products fail to exhibit the dominance of the divergent component in the stratosphere. Overall, stratospheric ERA5 and CFSv2 outperform the other reanalysis datasets, and only these two can reproduce the feature of the canonical kinetic energy spectrum with a distinct shift from a steeper slope (approximately −3) at lower wave numbers to a shallower slope (approximately −5/3) at higher wave numbers. In addition, the relative disparities among datasets increase dramatically with height, and they are more pronounced in the divergent component. It is also found that the correlations among these datasets are much weaker in the Tropics.

Evaluation of Multiweek Tropical Cyclone Forecasts in the Philippines

Abstract In the pursuit of providing tropical cyclone (TC) forecasts beyond the conventional time scales covered by weather forecasting in the Philippines, this study has examined the multiweek (i.e., from week 1 to week 4) TC forecast skill in the country. TC forecasts derived from three ensemble models, namely, the NCEP Climate Forecast System version 2 (CFSv2), the European Centre for Medium-Range Weather Forecasts Ensemble Prediction System (ECMWF), and the NCEP Global Ensemble Forecast System version 12 (GEFSv12) from 6 October 2020 to 31 October 2021 were verified. Results revealed that the ECMWF model is consistently the most skillful in multiweek TC prediction over the domain bounded by 110°–155°E and 0°–27°N in the western North Pacific. The ECMWF obtained hit rates ranging from 0.25 to 0.31, low false alarm rates of 0–0.33, and the highest equitable threat scores among the models. In contrast to this, the GEFSv12 and CFSv2 models had varying skills, with the former performing better in the first two weeks and the latter in longer lead times. It is further revealed that the three models generally underestimate the observed number of storms, storm days, and accumulated cyclone energy. Moreover, the study shows that the forecast TC tracks have a significant (p < 0.05) positional bias toward the right of observed tracks beyond week 1, and that they tend to propagate slower than observations especially in week 1 and week 2. These findings contribute to better understanding the strengths and limitations of these ensemble models useful for eventual provision of multiweek TC forecasts in the Philippines.

A procedure for obtaining forecast turbulence parameters in the atmospheric boundary layer for acoustic propagation

Propagation through turbulence causes a mean reduction in loudness of sonic booms and in some cases can cause significant fluctuations in loudness about that mean. Understanding these effects is important for accurately planning loudness levels during upcoming community noise tests with the X-59 aircraft as well as mitigating the risk for excessive loudness due to turbulence. Current methods for modeling acoustic propagation through turbulence require several atmospheric parameters including the atmospheric boundary layer (ABL) height, friction velocity, mixed-layer velocity scale, surface-layer temperature scale, as well as ambient pressure, temperature, and humidity within the ABL. Accurate forecasts of these parameters are needed for X-59 flight planning and may also be useful for planning future supersonic aircraft certification flight tests. This presentation showcases one method for obtaining the forecast turbulence and ambient atmospheric parameters from freely available Climate Forecast System Version 2 data. Forecast parameters are compared to measurements to assess accuracy and utility.

Simulations of X-59 sonic thumps and traditional sonic booms propagated around the world for three atmospheric models

Propagation simulations of sonic booms from supersonic aircraft through atmospheric data over time at fixed locations provide the opportunity to assess noise exposure statistics for different climate regions. Knowledge of climate-based differences in sonic boom noise exposure statistics is important to ensure that future civil supersonic aircraft noise certification standards are globally applicable and effective. In this presentation, simulated sonic booms from the NASA X-59 Quesst quiet supersonic aircraft and conventional supersonic aircraft were propagated through atmospheric data at 100 locations across the world using PCBoom. Noise exposure statistics are compared for propagation results from three different atmospheric databases (NOAA Global Forecast System, NOAA Climate Forecast System Version 2, and the ECMWF Reanalysis Version 5). These atmospheric models were chosen due to their global coverage, popularity, and database availability. Preliminary statistical models are fit to assess the impact of several factors including flight direction, season, ground elevation, and climate on noise exposure size and loudness. Areas with prevalence of higher noise due to their climate are identified, which could help inform future supersonic aircraft noise standards.

Statistical Evaluation of the Performance of Gridded Daily Precipitation Products from Reanalysis Data, Satellite Estimates, and Merged Analyses over Global Land

The Beijing Climate Center of the China Meteorological Administration (BCC/CMA) has developed a gauge-satellite-model merged gridded daily precipitation dataset with complete global coverage, called BCC Merged Estimation of Precipitation (BMEP). Using the unified rain gauge dataset from the CPC (CPC-U) as the independent benchmark, BMEP and the four most widely used global daily precipitation products, including the Global Precipitation Climatology Project one-degree daily (GPCP-1DD), the NCEP Climate Forecast System Reanalysis (CFSR), the Interim ECMWF Re-analysis (ERA-interim), and the 55 year Japanese Reanalysis Project (JRA-55), are evaluated over the global land area from January 2003 to December 2016. The results show that all gridded datasets capture the overall spatiotemporal variation of global daily precipitation. All gridded datasets can basically capture the overall spatiotemporal variation of global daily precipitation. However, CFSR data tend to overestimate precipitation intensity and exhibit a spurious positive trend after 2010, attributed to the transition from CFSR to NCEP’s Climate Forecast System Version 2 (CFSv2). On the other hand, JRA-55 and ERA-interim data demonstrate higher skill in characterizing spatial and temporal variations, bias, correlation, and RMSE. GPCP-1DD data perform well in terms of bias but show limitations in detecting the interannual variability and RMSE of daily precipitation. Among these evaluated products, BMEP data exhibit the best agreement with CPC-U data in terms of the spatiotemporal variation, pattern, magnitude of variability, and occurrence of rainfall events across different thresholds. These findings indicate that BMEP gridded precipitation data effectively capture the actual characteristics of daily precipitation over global land areas.

Improving the CFSv2 Seasonal Precipitation Forecasts across the United States by Combining Weather Regimes and Gaussian Mixture Models

Abstract Although seasonal climate forecasts have major socioeconomic impacts, current forecast products, especially those for precipitation, are not yet reliable for forecasters and decision-makers. Here we developed a novel statistical–dynamical hybrid model for precipitation by applying weather regimes (WRs) and Gaussian mixture models (WR-GMM) to the National Oceanic and Atmospheric Administration’s Climate Forecast System, version 2 (CFSv2), precipitation forecasts across the continental United States. Instead of directly forecasting precipitation, WR-GMM uses observed precipitation from synoptic patterns similar to the future CFSv2 forecast. Traditionally K-means has been used to classify daily synoptic patterns into individual WRs, but the new GMM approach allows multiple WRs to be represented for the same day. The novel WR-GMM forecast model is trained on daily Climate Forecast System Reanalysis (CFSR) geopotential height and observed precipitation data during the 1981–2010 period and is verified for 2011–22. Overall, the WR-GMM method outperforms the CFSv2 ensemble forecast precipitation in terms of root-mean-square error and for Pearson correlation coefficient for lead months 1–4. Previous studies have used global climate models to forecast WRs in the Pacific Ocean and Mediterranean Sea regions, usually with an emphasis on winter months, but the WR-GMM model is the first of its kind that promises great untapped potential to improve precipitation forecasts produced by CFSv2 across the continental United States.

Review article: Towards strongly coupled ensemble data assimilation with additional improvements from machine learning

Abstract. We assessed different coupled data assimilation strategies with a hierarchy of coupled models, ranging from a simple coupled Lorenz model to the state-of-the-art coupled general circulation model CFSv2 (Climate Forecast System version 2). With the coupled Lorenz model, we assessed the analysis accuracy by strongly coupled ensemble Kalman filter (EnKF) and 4D-Variational (4D-Var) methods with varying assimilation window lengths. The analysis accuracy of the strongly coupled EnKF with a short assimilation window is comparable to that of 4D-Var with a long assimilation window. For 4D-Var, the strongly coupled approach with the coupled model produces more accurate ocean analysis than the Estimating the Circulation and Climate of the Ocean (ECCO)-like approach using the uncoupled ocean model. Experiments with the coupled quasi-geostrophic model conclude that the strongly coupled approach outperforms the weakly coupled and uncoupled approaches for both the full-rank EnKF and 4D-Var, with the strongly coupled EnKF and 4D-Var showing a similar level of accuracy higher than other coupled data assimilation approaches such as outer-loop coupling. A strongly coupled EnKF software framework is developed and applied to the intermediate-complexity coupled model SPEEDY-NEMO and the state-of-the-art operational coupled model CFSv2. Experiments assimilating synthetic or real atmospheric observations into the ocean through strongly coupled EnKF show that the strongly coupled approach improves the analysis of the atmosphere and upper ocean but degrades observation fits in the deep ocean, probably due to the unreliable error correlation estimated by a small ensemble. The correlation-cutoff method is developed to reduce the unreliable error correlations between physically irrelevant model states and observations. Experiments with the coupled Lorenz model demonstrate that strongly coupled EnKF informed by the correlation-cutoff method produces more accurate coupled analyses than the weakly coupled and plain strongly coupled EnKF regardless of the ensemble size. To extend the correlation-cutoff method to operational coupled models, a neural network approach is proposed to systematically acquire the observation localization functions for all pairs between the model state and observation types. The following strongly coupled EnKF experiments with an intermediate-complexity coupled model show promising results with this method.

Prediction of Intermonth Modes of Winter Air Temperature over China

Abstract The features and causes of the leading intermonth modes of winter surface air temperature anomalies (SATA) over China are investigated, and associated prediction models are developed. The first three intermonth modes of winter SATA over China are obtained by extended empirical orthogonal function (i.e., EEOF1–3) analysis. The results show that EEOF1 represents consistent variations in the whole winter, with a variance contribution of 32.3%, whereas EEOF2 and EEOF3 show spatiotemporally inconsistent changes, and their variance contributions are 16.9% and 12.5%, respectively. EEOF2 has out-of-phase variations between December and January–February, and EEOF3 exhibits a temporal warm–cold alternating pattern, with spatially reversing changes over northwestern and southern China. However, the Climate Forecast System, version 2 (CFSv2) presents a limited prediction skill for winter SATA over China and their intermonth modes. Further investigations indicate that the September sea ice over the Barents–Laptev Seas, the November snow cover over western Europe and East Asia, and the November northern Atlantic sea surface temperature can be, respectively, adopted to develop prediction schemes for the consistent mode (EEOF1 scheme) and two inconsistent modes (EEOF2 and EEOF3 schemes) based on specific mechanisms. These schemes show effective performances in predicting both individual modes and the reconstruction field of SATA over China. The temporal correlation coefficients (TCCs) between cross-validation results and observations are 0.48, 0.51, and 0.31 for the EEOF1–3 modes, respectively (the 90% confidence level is 0.27). For the reconstruction field, the TCCs are 0.40, 0.27, and 0.45 in December, January, and February, respectively, which are much higher than those of the CFSv2 outputs (0.23, −0.16, and −0.09).

Abrupt sub-seasonal rainfall variability over India during summer monsoon 2021: Interaction between midlatitude and tropical circulation

The Indian subcontinent witnessed abnormal subseasonal rainfall variations during the summer monsoon season 2021 with abrupt change (54%) in rainfall from August to September. According to India Meteorological Department, in August Indian landmass received 81% rainfall of its long period average (LPA), which is lowest since 2002. On the other hand, anomalously high rainfall (+ 35% of its LPA) was received in September. The physical mechanisms responsible for such unusual changes in monthly/subseasonal rainfall are unravelled in this study. In August 2021, the southward displacement of Asian jet (SWDAJ) corroborated by westward and northwestward extension of the Western north Pacific (WNP) anticyclone induced moisture divergence over the monsoon trough region and northwestern parts of India are mainly responsible for the deficit rainfall. In the East Asian region, strong low-level westerly wind anomalies in the northern flank of WNP anticyclone couples with the upper-level westerlies associated with SWDAJ and the resulting interaction due to the barotropic–baroclinic coupling maintains WNP circulation. Apart from WNP anticyclone, anomalous Sea Surface Temperature warming over the southeast equatorial Indian Ocean also contributed to the suppressed rainfall over India. In contrast to August, during September the northward displacement of the Asian jet and associated anomalous upper-level anticyclone over northwest India linked to the mid-latitude wave train emanating from the northeast Atlantic contributed to enhanced rainfall over northwest India. Anomalous upper-level high over central Asia/northwest of India increased the moist baroclinic instability over northwest India, and enhanced rainfall there. The Climate Forecast System version 2 (CFSv2) model based on May initial conditions is some extent predicted rainfall change from August and September 2021 over India with much lesser magnitude. Inability of the model in capturing the changes in Asian jet, midlatitude wave train and associated teleconnections/interactions to Indian Summer Monsoon (ISM) rainfall limits its prediction skill. Further, the North American Multi Model Ensemble (NMME) models also could predict the correct sign of ISM rainfall anomaly in August and September, with May initial conditions, but failed to capture the spatial patterns, which needs further investigation.